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Abstract 



Using a specific model for tlie expansion rate of the Universe as a function 
of scale factor, it is demonstrated that the equation of state of the dark 
energy cannot be determined uniquely from observations at redshifts z ^ 
a few unless the fraction of the mass density of the Universe in nonrelativistic 
particles, i^M, somehow can be found independently. A phenomenological 
model is employed to discuss the utility of additional constraints from the 
formation of large scale structure and the positions of CMB peaks in breaking 
the degeneracy among models for the dark energy. 
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I. INTRODUCTION 



Although observations of acoustic peaks in the Cosmic Microwave Background (CMB) 
fluctuation spectra are sensitive to the curvature of spacetime, and appear to require some 
form of cosmological dark energy they may be less useful for discriminating among 

different equations of state for the dark energy. Various studies (e.g. have concluded 

that the signatures of different types of spatially smooth, evolving dark energy on the CMB 
spectra are relatively undiscriminating, in part because CMB fluctuations were formed at 
high redshift, when the dark energy was not a prominent constituent, although spatial 
fluctuations in the dark energy could lift the degeneracy considerably 0. By contrast, 
observations of sources at small to moderate redshifts (2; ^ a few) would probe epochs 
where the dark energy is prominent, leading to suggestions that observations of Type la 
supernovae pO|JTl|JT3H T7| or of galaxy counts |T9|-pT| could be used to determine the nature 
of the dark energy most effectively. However, the efficacy of these observational programs 
is controversial, as there is considerable degeneracy among the predictions of different dark 
energy models for, for example, the luminosity distance as a function of redshift ||18|| . 

One approach to the analysis of low redshift data would be to presume particular classes 
of models, perhaps parametrized by a series expansion of w{z) = p{z) / p{z), and then attempt 
to constrain the model parameters by a likelihood or Bayesian method (e.g. [^^-|T^,|TB|]). A 
second approach is to fit the data by a parametric representation of, for example, the lu- 
minosity distance as a function of redshift, and then analyze the results to constrain the 
properties of the dark energy indirectly (e.g. p!3|-p!7|). Although the first approach is prefer- 
able (because it generally avoids difficulties associated with differentiating data, and allows a 
simpler assessment of uncertainties in derived parameters, in addition to outlining its implicit 
assumptions more clearly), the second approach is more useful for understanding whether 
or not such observations can ever yield a unique solution for the dark energy equation of 
state. 

In this short note, the question considered is: Suppose one could analyze data from 
z ^ a few to determine the expansion rate of the Universe as a function of scale factor, 
H[a), exactly. Could one then determine the equation of state of the dark energy exactly 
from these measurements alone? The answer, as we shall see, is no: a separate determination 



of the fraction of the Universe in the form of nonrelativistic particles is needed |22 



II. A SIMPLE ILLUSTRATIVE MODEL 

The necessity of separately determining the density of nonrelativistic particles in order 
to deduce the equation of state of the dark energy can be demonstrated most easily by 
constructing a specific example. Suppose the data show that the Universe is expanding 
according to the simple relationship 

^ HI (n. ^ ^) . (1) 

where a = (1 + ^)^^. One (tempting) interpretation of this result would be that the Universe 
consists of two components, nonrelativistic matter contributing a fraction of the closure 
density, and a cosmological constant contributing a fraction f2i = 1 — ^2- However, this 
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interpretation is not unique: there are models involving a scalar field with nonconstant 
effective potentials ^(0) that can lead to Eq. (0). 

To see this, we can construct an explicit model with both a scalar field and nonrelativistic 
matter, so that 

with Qm 7^ Q2 in general. The scalar field energy density is therefore 

where po = SHq/SttG. Differentiate with respect to a to find 

3po{n2 - ^Im) _ 3(p0 + P^) 



da a 
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(4) 



where the second equality follows from conservation of energy for the scalar field. Thus, we 
find simply that = —poQi, and 

^' = Ps + Ps = ^^^^ " ^''^ 



V{<^) = ^— = PO [n^ + j • (5) 

Note that the first of Eqs. (|^) requires that Q2 > ^m- We can determine V{(f)) by solving 
the first of Eqs. (^) for 0(a), inverting to find a(0), and substituting into the second of 
Eqs. (H). Combining Eqs. and the first of Eqs. implies 



Define 



then if ip = 



^= aiid = ; (7) 



1 

V c/a; / 1 + x ^ ^ 



which has the general solution 



^ = 7/'oo±2sinh-i f^j , (9) 



where -000 is a constant of integration. The sign in this solution can be chosen arbitrarily; 
chosing the — sign gives a solution in which -0 increases with increasing a. If we let ip = i/jq 
today, then we find 
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ip = ipo + sinh 
which we may invert to find 



1^' 2sinh-i I 



(10) 



1 



cosh f — I —J— sinh f — 



■^ 2 



(11) 



where i/j = ijj — i/jq. The second of Eqs. (||) then determines the potential to be 

n 2 



cosh f — I —J— sinh f — 



(12) 



Note that -ip —>■ 2 coth ^ y 1/^72 and 1^(0) ^ Po^i asymptotically. Also, when ^2 = 1 — ^1"^ 

1, the potential becomes exponential: V{(j)) — > |po(l — ^M)e~'^. 

The key question here revolves around prior assumptions. The "simplest" interpretation 
of Eq. (|1]) is that the energy density of the Universe is the sum of contributions due to 
nonrelativistic matter and a cosmological constant. However natural this interpetation may 
seem to be, it only follows if one assumes either that the dark energy is in the form of a 
cosmological constant, or that the mass density in nonrelativistic form has closure parameter 
^2- If one is somewhat more agnostic on either of these points, then a range of possibile 
interpretations accounts for the data equally well. The degeneracy is not just a property of 
models for which H^{a) is given by Eq. (|1]). If the data were to favor a model of the form 



H\a) = Hi 



(13) 



with g{a) some function, and VLig{l) + ^^2 = 1? then without additional prior assumptions 
the most one can say for sure is that = po [flig{a) + {fl2 — ^Af)/ct^], which leaves open 
the possibility of a range of different solutions for V{(j)). 

This model illustrates that if one makes no prior assumptions about the form of the 
quintessence potential V{(j)), then it is impossible to fit low redshift data alone to find 
V{(f)) uniquely unless Qm is known independently - in this example, a range of models, 
parametrized by the value of / = will fit equally well. Constraints on the value of / 

limit the range of acceptable models. Although minimally / < 1, we know that this model 
would be unable to account for the formation of large scale structure if / is too small. If 
we assume that only the nonrelativistic matter can clump (i.e. that the dark energy and 
relativistic matter remain smooth), then, in linear theory, the density contrast D{a) obeys 
the equation (e.g. p^) 




4:nGpM 



D = 



(14) 



Eq. (|T^ has both growing and shrinking modes, -D+(a) and D-{a), respectively. If we 
include radiation (which is a known component of the Universe, and subdominant today, 
but important in determining linear theory growth rates), then Eq. (|l|) must be altered to 
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where Q1 + Q2 + ^rad = 1 and, for three low mass or massless neutrinos, Qradh^ 7 ~ 8.5 x 10~^ 
if ho J = Hq/{70 km s^^ Mpc"^). The growing mode mainly amplifies perturbations between 
tteq = VLrad/^2 and tti = (^2/^1)^^^- During this phase, oc a~^, and -D+(a) oc a"+ , where 

^+(/) = ^(\/l + 24/-l). (16) 

The overall linear theory growth factor is ~ {ai/ aeqY+'^^\ so the ratio of the growth factor 
for / 7^ 1 to its value for / = 1 is approximately 



where the approximation assumes 1-/^1. Thus, we estimate that for the growth factor 
to be within a factor of two of its value for / = 1, we must have 1 — / < 51n2/31n(ai/aeg) — 
0.12. 

The solid line in Fig. 1 shows the ratio of the growth factor for VLm < ^2 to its value for 
0.2 = ^A/, that is 

as a function of Om for O2 = 0.5, where i5+(a; ^a/, ^2) is the linear theory growth factor 
evaluated at scale factor a, given Om and O2. These results were found by solving Eq. 
(p!^) numerically from a -C agq to the present day, assuming identical initial perturbation 
amplitudes. (Relativistic particles were included as a smooth component via Eq. \T^\.) The 
numerical results show that for Rsp/p > 0.5, then we must require Om ~ 0.44, corresponding 



to / ^ 0.88, approximately the same as the analytic estimate of the previous paragraph. 
(The estimate is accurate to about 20% for O2 ^ 0.1.) 

Linear theory also predicts velocity fluctuations proportional to the peculiar velocity 
factor T^{a; Om, O2) = dlnD^{a; O^j, 02)/dlna. The dashed line in Fig. 1 shows 

the ratio of the peculiar velocity factor at a = 1 for Om < ^2 to its value for Om = O2. 
Although there is some dependence of Rgv on ^2^,/, it is not as extreme as for Rsp/p- 

Thus, primarily from the requirement that linear density perturbations can grow sub- 
stantially from small values in the early Universe, we know that if Eq. ([^) were truly exact, 
/ must be close to one. Imposing limits on the range of plausible Om from independent, 
physical and phenomenological arguments amounts to using prior information to restrict 
the possibilities to be compared with the data from observations at z ^ a few. (The use 
of prior information from large scale structure and CMB anisotropies has been discussed in 
e.g. [p!5| , p5| -[27|| . ) This does not eliminate the degeneracy implicit in Eq. (||) altogether, but 



does diminish its significance by constraining the range of acceptable values of / severely. 
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However, the fact remains that if one were to try to determine the equation of state of the 
dark energy from observations that, hypothetically, yield Eq. (|l]), then one can never hope 
to find a unique result from an analysis that ignores constraints from the development of 
large scale structure (or any other considerations that yield independent information about 
Qm, such as peculiar velocities, CMB fluctuations and weak lensing). Moreover, even if 
such restrictions are imposed, unless Qm is determined precisely by them, some degeneracy 
remains in the construction of V{(f)) from observations at few. 




FIG. 1 Relative growth factor (solid) and velocity factor (dashed) 
of linear perturbations for models with Qm < Q2 = 0.5 that expand 
according to Eq. (|l|). 



III. CAN Qm and n2 DIFFER SUBSTANTIALLY? 

If Eq. (|l]) with ^2 substantially different from Q^j is an excellent approximation at low to 
moderate redshifts, it can never he exact. In this case, there must be a transition in H'^{a), so 
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that the component of the density that scales hke is predominantly due to nonrelativistic 
particles for a period sufficient to grow large scale structure. The implicit deviations from 
Eq. (|I|) would contain information on VLm as well as f22, and probably additional parameters 
as well. In this situation, it might be possible to ascertain the equation of state of the dark 
energy from measurements at low to moderate redshifts alone, provided that the imprint 
of VLm and the additional parameters in H'^{a) can be discerned from these observations. 
The question is, if Eq. (|I|) is an accurate (but inexact) fit at 2; ^ a few, and fl2 differs 
substantially from Qm, how well could we determine ^Im, and hence break the degeneracy, 
observationally? 

To make these issues more concrete, suppose in reality Eq. (|I|) is actually the large a 
limit of a more general, exact relationship, say 



H\a) = Hl{a)F [ — )+ Hl{a) 



at 



(20) 



where F[q) — as g ^ and F{q) ^ 1 as g ^ 00, with < 1 a transition value of the 
scale factor, H^{a) given by Eq. (P, and, for example, 

Hl{a) = H^,(n, + ^) ; (21) 



since the model is fiat, 



i = {n^ + n2)Fo + {n< + nM){i - i^o), (22) 



where Fq = F{l/at). Such a model could be consistent with the growth of large scale 
structure, but might still be hard to distinguish from Eq. (|l|) based on observations at 
low to moderate redshift, depending on the value of at, and the form of F{a/at). Eq. 
(^), although admittedly ad hoc, is useful for studying the extent to which observations of 
various sorts could discern that Q2 7^ ^m, and pin down the difference sufficiently to break 
the degeneracy discussed above. 

For a fiat cosmology, Eq. (^0]) implies that the radial coordinate for a source at a = 



l + z)-' is 



rs{a; Qm, ^2, F) = //(I) - r]{a) 



1 db 



a b^Hib) 



^-1 dh 

° -^^ Jb{nj)^ + 1 - fie + [(^^< - ^iW + - ^2][F^ - F{b/at)]} 



where fig = ^1-^0 + fi<(l ~ -^0)5 and dr] = dt/a = da/a^H{a). (Here, the notation 
rs{ci',^M)^2-, F) means the radial coordinate as a function of a given VLm and ^2 > ^m, 
as well as a set of parameters required to define F{a/at).) If F{a/at) does not vary much 
over the range of a covered by the observations, then the data would only determine Qg 
accurately, and would yield little useful information on ^2 or Qm'- only the combination 
fl2Fo + fiA/(l — Fq) = 1 — Qe could bc deduced. This situation is possible provided that at 
is fairly small, below the range covered by observations, and F{a/at) varies relatively slowly 
as a function of a/at for large values of a/at- In this case, we cannot expect to eliminate 



7 



the degeneracy among models directly from observations of low to moderate redshift sources 
alone. Constraints from other sorts of observations that probe Qm directly would be needed 
to deduce the quintessence model as accurately as possible. 

More quantitative statements only can be made in terms of specific F{a/at), and require 
some exploration of parameters. Let us examine one choice: fii + ^2 = 1 = f2< + ^m, in 
which case Eq. ([2^ ) is satisfied exactly, and Eq. (p3D becomes 

rs{a;n2,nM,F)=H,' (24) 

^b{n, + (1 - n^)^ - (^2 - ^^m)(1 - b^) [1 - F{b/at)]} 




z 



FIG. 2 Magnitude differences among models with = 0.5 for 
at = 0.2 (solid) and at = 0.5 (dashed) as functions of z for, from 
the top for each at, flu = 0.1, 0.2, 0.3 and 0.4. 

Eq. ( pH ) may be used to compute magnitude differences among models with given Q2, 
VIm and F(a/at); for illustrative purposes, let us choose 1 — F{a/at) = e~"/"'. Fig. 2 depicts 
the magnitude difference 
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Am{z;il2:^M:at) = 5.0 log 



10 



rs{{l + z) ^■,nM,^2,at) 
rs{{l + z)-^;n2,n2,at) 



(25) 



relative to the model with = = 0.5 as a function of z assuming at = 0.2 (solid) and 
at — 0.5 (dashed) for (from top to bottom in each case) — 0.1,0.2,0.3 and 0.4. Not 

surprisingly, the magnitude differences are smaller for the larger value of at, and decrease as 
Qm ^2- Nevertheless, the largest magnitude difference shown for 2; < 3 is 0.29 mag, for 
{z, at, flu) = (3, 0.5, 0.1); for {z, at, fin) = (3, 0.2, 0.1), the difference is 0.077 mag. 




FIG. 3 Fractional deviations of the angular size of the (sonic) 
horizon at a = 10~^ from the = ^2 — 0.5 model as functions 
of Qm < ^^2- The sohd hue is for at — 0.2, and the dashed hue for 
at = 0.5. For computing the angular sizes. Ho — 70 km s~^ Mpc~^ 
was adopted. 



Fig. 3 shows results for 
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r/(10 ^;l]A/,1^2,Qt)r5(10 ^1^2,^2,af) 



-1 , 



(26) 



the fractional difference between the angular size of the (sonic) horizon at a = 10"^ as a 
function of Vt^ < ^2 and its value for VLm = ^2 for the two different values of at. (Deviations 
of the sound speed from 1 / have not been included.) Larger fractional deviations arise for 
smaller VLm and a^. The largest deviation shown in the two figures is 12.7% for (^2^,04) = 
(0.1, 0.2); the largest fractional deviation for at = 0.5 is about 1.3%, and occurs near Qm = 



Although this is only one example, based on a particular {ad hoc) choice of F{a/at), 
it illustrates that while data at both low to moderate z as well as CMB observations can 
yield information that lifts the degeneracy between Q2 and ^Im, considerable precision may 
be needed to determine the equation of state of the dark energy accurately. Note that no 
attempt has been made to fit simulated data based on Eq. (p6|) to figure out how well H{a) 
could be determined in practice; the comparisons made in Figs. 2 and 3 are based on exact 
functions. The observational challenge of distinguishing among models at these levels could 
be even greater than these relatively small differences would indicate. On the other hand, 
Figs. 2 and 3 were constructed based on a particular {ad hoc) model for H'^{a) designed for a 
phenomenological study of how easily information on VLm could be gleaned from luminosity 
distances and CMB acoustic peaks. The task could be easier or harder in cosmologies based 
on particular physical theories for the dark energy, depending on the details of the model. 



To conclude, the point of this paper is not to suggest that it is absolutely inconceivable 
that the equation of state of the dark energy could be measured, but rather to show, via 
a specific example, how measurements at 2; ^ 3 admit a continuum of interpetations in 
terms of evolving dark energy fields unless Qm can be determined separately somehow. The 
specific example, Eq. (|I|), can be intepreted "naturally" in terms of a Universe containing 
a mixture of nonrelativistic particle dark matter plus a cosmological constant, but can be 
interpreted equally well in terms of a quintessence model with a range of potentials, Eq. (|T2|). 
This specific example sheds some light on why simulated analyses of low to medium redshift 
data appear to allow degenerate interpretations: it is always possible that the component 
of the total energy density of the Universe that evolves like a~^ is only partially due to 
nonrelativistic particles, with the rest arising from quintessence. Unique interpretations 
can only be obtained in terms of specific models for the quintessence that forbid such a 
conspiracy a priori, by either specifying the form of V{(f)) or fixing the value of Qm- 

Realistically, the extent to which the p oc a~^ constituent must be due to nonrelativis- 
tic particles is constrained by the requirement that large scale structure formation evolve 
"normally," that is, unimpeded by the existence of a component with density proportional 
to that is incapable of clustering. If Eq. (|I|) were truly exact, then fl2 would have to 
be close to Qm for large scale structure to grow, but even in this case, a limited degeneracy 
remains in the determination of V{(j)). If ^2 and Qm differ substantially, then Eq. (|T|) cannot 
be exact, raising the question of how well one could discern the two parameters separately 
from observations. To examine this issue, a modified (phenomenological) expansion law. 



0.225. 



IV. DISCUSSION 
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Eq. (pOD, that encodes information on both ^2 and Qm was introduced. In the context of 
a particular model that is consistent with large scale structure formation, but reduces to 
Eq. (|^) with Q2 > at large enough a, we have seen that information about Qm could 
be gleaned both from measurements of the positions of CMB peaks and from luminosity 
distance determinations. However, the deviations among models with various Qm < ^2 
can be quite small, which would still pose a substantial challenge for programs that aim to 
determine the equation of state of the dark energy. 

Of course, in the context of specific models for the quintessence field, embodied in partic- 
ular forms for V{(j)), the analyses may not encounter pronounced degeneracies. At present, 
though, there is little compelling reason to assume any particular V{(j)) a priori. When 
analyzing data for a particular form of V{(j)), one is primarily engaged in estimating the pa- 
rameters of the model, as well as Qm- Thus, one might be able to find the best fit model of 
a particular type (e.g. V{(j)) = constant), without being able to tell if the underlying model 
would be favored by the data if other possibilities were admitted. If the allowed ranges 
of parameters for a particular quintessence model do not shrink with the accumulation of 
data, one can conclude that the model is inadequate with confidence. However, there is no 
guarantee that a given model is correct even if the data seem to converge on a unique set of 
parameters, as Eqs. (|ip and (|r^) illustrate. 
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